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The decays of J/ip and ip{2S) to nK'gK + c.c. are observed and measured for the first time, and the perturbative 



1. Introduction 

Since the discovery of the J/il> at 
Brookhaven [1] and SLAC [2] in 1974, more 
than one hundred exclusive decay modes of the 
J/ij) have been reported. According to Ref. [3], 
direct hadronic, electromagnetic and radiative 
decays make up roughly 65%, 14%, and 7% of 
the total J/ij.) decay width, respectively. How- 
ever, the measured hadronic decay channels sum 
up to less than 35 %. The BESII data sample of 
5.8 xlO^ J /ijj events provides a good opportunity 
to search for missing J /ijj hadronic decays. 

In 2004, BESII reported the observation of 
an enhancement X(2075) near the threshold of 
the invariant mass spectrum of pA in J /tp -—^ 
pK~A decays. The mass, width, and product 
branching fraction of this enhancement are M = 
2075 ± 12 (stat.) ± 5 (syst.) MeV/c^, T = 90 ± 
35 (stat.) ±9 (syst.) MeV/c^ [4], and B{J/^Jj 
K-X)B{X ^ pA+c.c.) = (5.9±1.4±2.0)xl0-^ 
respectively. The study of the isospin conjugate 
channel J/ip nKgA is therefore important not 
only in exploring new decay modes of J/tp but 
also in understanding the X(2075). 

The 5.8 X 10^ J/^p and 1.4 x 10^ ^'(25') events at 
BESII also offer a unique opportunity to search 
for new decay modes of J/ip and tp{2S) and test 
the "12% rule" in hadronic decays. In perturba- 
tive QCD, hadronic decays of the J/tp and tp{2S) 
are expected to proceed dominantly via three glu- 
ons or a single direct photon with widths propor- 
tional to the square of the cc wave function at the 
origin, which is well determined from dilepton de- 
cays. Thus for any hadronic final state h, the J/tp 
and tp{2S) decay branching fractions should sat- 
isfy the so called "12% rule" [5]. 

_ B{i;{2S) ^ h) ^ B{i,{2S) ^ e+e-) ^ 
B{J/iP^h) ~ B{J/iP^e+e-) ~ ' 

The leptonic branching fractions are taken from 
the particle data group (PDG) [6] tables. It 
is roughly obeyed for a number of exclusive 
hadronic decay channels except some FP, PP 
and VT channels [7] [8] [9], where P, V and T de- 
note members of the pseudoscalar, vector and 
tensor nonets, respectively. 

In this paper, the first observation and mea- 



surement of J/V' and '(p{2S) to nK^A. -\- c.c, as 
well as a test of the perturbative QCD 12% rule 
are presented. The A* and N* resonance struc- 
tures in J/tp ^ nKgA+c.c. are also shown, where 
no obvious enhancement near nA threshold is ob- 
served. The upper limit on the branching fraction 
of J/tp ^ ^s^' ^ ~^ determined. 



2. The BES Detector 

The upgraded Beijing Spectrometer detec- 
tor (BESII) is located at the Beijing Electron- 
Positron Collider (BEPC). BESII is a large solid- 
angle magnetic spectrometer which is described 
in detail in Ref. [10]. The momentum of charged 
particles is determined by a 40-layer cylindri- 
cal main drift chamber (MDC) which has a mo- 
mentum resolution of (jp/p=1.78%-\/T+p^ jjj 
GeV/c). Particle identification is accomplished 
using specific ionization (dE/dx) measurements 
in the drift chamber and time-of-flight (TOE) in- 
formation in a barrel-like array of 48 scintilla- 
tion counters. The dE/dx resolution is OdE/dx — 
8.0%; the TOE resolution for Bhabha events is 
ctof = 180 ps. Radially outside of the time- 
of-flight counters is a 12-radiation- length barrel 
shower counter (BSC) comprised of gas tubes in- 
terleaved with lead sheets. The BSC measures 
the energy and direction of photons with resolu- 
tions of (je/E ~ 21%/B {E in GeV), = 7.9 
mrad, and Cz = 2.3 cm. The iron flux return 
of the magnet is instrumented with three double 
layers of proportional counters that are used to 
identify muons. 

A GEANT3 based Monte Carlo (MC) program 
(SIMBES) [11] with detailed consideration of the 
detector performance is used. The consistency 
between data and MC has been carefully checked 
in many high purity physics channels, and the 
agreement is reasonable. More details on this 
comparison can be found in Ref. [11]. The detec- 
tion efficiency and mass resolution for each decay 
mode in this analysis are obtained from MC sim- 
ulation. 
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3. Analysis 



The analyzed J/iJj and xp{2S) nK'^A with 
r+TT" and A pn^ (and c.c.) final states 
contain four charged tracks and an undetected 
neutron or anti-neutron. We require the candi- 
date events to satisfy the following common se- 
lection criteria: 

1. Events must have four good charged tracks 
with zero net charge. A good charged track 
is a track that is well fitted to a three- 
dimensional helix, originates from the inter- 
action region and has a polar angle 9 in the 
range | cos9\ < 0.8. Because of the long de- 
cay lengths before Kg and A decay to tt+tt" 
and p7r+, the interaction region is defined 
as Rxy < 0.12 m and \z\ < 0.3 m. Here, 
Rxy is the distance from the beamline to 
the point of closest approach of the track to 
the beamline, and \z\ is the distance along 
the beamline to this point from the interac- 
tion point. 

2. For each charged track in an event, x^pioi^) 
is determined using both dE/dx and TOF 
information: 

Xpioii) = xls/dxii) +XTOFii)^ 

where i corresponds to the particle hypoth- 
esis. A charged track is identified as a p if 
XpiD fo'" tti^ P hypothesis is less than those 
for the 77 or if hypotheses. For the channels 
studied, one charged track must be identi- 
fied as a p or p. 



3.1. Measurement of J ftp — »■ nKgA + c.c. 

For J/tp nK^A -> pmr+-n:-7T+ , the ^ 
TT+TT" and A — > p7r+ decays are reconstructed us- 
ing secondary vertex fitting, and the 7r+ from the 
A decay is identified. To select A(A) [A(A) mass 
selection], |Mp„-(p^+) - 1.115| < 0.012 GeV/c^ is 
required, and to select Kg {K^ mass selection), 
\M„+„- - 0.497| < 0.02 GeV/c^ is required. To 
reject the backgrounds from channels containing 
a K^ but no A, like e.g. J/ijj pK'^Y,^ + c.c, 
we require Lxy{A), the distance from the recon- 
structed A vertex to the event origin, to be larger 
than 5 mm. 



Figure 1 is the missing mass spectrum de- 
termined from the charged tracks in J/ tp — > 
TiKgA + cc. candidate events satisfying A(A) and 
Kg mass selections and Lxy{A) > 5 mm. A clear 
peak at the nominal neutron mass is observed. 
The second peak in the high missing mass re- 
gion comes from J/tp nKgY."{1385) + c.c. and 
J/ip S~E+(1385) -I- c.c. backgrounds. To sup- 
press background and improve the resolution, a 
one constraint (IC) kinematic fit with a missing 
neutron is applied under the J/ip pmT'^'jT~w~^ 
hypothesis. The distribution of IC fit x|„t+^-^+ 
for the above selection is shown in Fig. 2. The 
agreement between data and MC simulation is 
reasonable, and in the following, xic < 5 is re- 
quired. 
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Figure 1. The missing mass determined from the 

charged tracks in J/tp nK^A + c.c. candidate 
events satisfying A (A) and Kg mass selections 
and Lxy{A) > 5 mm. 



Figure 3 shows the scatter plot of TOp„+ ver- 
sus m^+^- , and clear A and Kg signals are seen. 
Figure 4 shows the A decay length distributions 
for data and MC simulation for events satisfy- 
ing A (A) and Kg mass selection requirements 
and Xic < ^- '^^^ missing mass distribution 
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Figure 2. The distributions for the IC fits 
to the J pmr~^ Tr~ tt'^ hypothesis for J/ip —^ 
nKgA candidate events for the selection used in 
Fig f . The crosses are data; the full histograms 
are the sum of MC simulation of J/t/j nKgA 
and background determined from Kg sidebands 
(0.06 < |M(7r+7r-) - 0.497| < 0.08 GeV/c^). 



of charged tracks for events satisfying these re- 
quirments plus L^yiA) > 5 mm is shown in Fig. 5, 
and a very clean neutron peak is seen. A fit with 
a Gaussian function yields a mass value consis- 
tent with that of the neutron. The tt+tt" invari- 
ant mass spectrum is shown in Fig. 6, and a Kg 
signal is clearly seen. 
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Figure 3. Scatter plot of mp^+ versus m^+j^- for 
J/tp ^ TiKgA candidate events satisfying Xic < 
5 and Lxy{A) > 5 mm. 



Figure 4. The A(A) decay length distributions 
with the A(A) and Kg mass selection require- 
ments and Xic < 5 for data and MC. The his- 
togram is the sum of signal MC and background 
from Kg sidebands, and the crosses are data. 



Backgrounds for J/tp nKgA which con- 
tribute to the peak in the Kg signal region mainly 
come from J/ijj nK^f]° and J/ip P-ftTgS" 
that survive selection criteria. Normalizing with 
the corresponding branching fractions and the 
number of J/tp events in the data sample, a to- 
tal of 42 ± 8 J/xp nK%fP and 12 ± 3 J/ip 
pK^gT," background events are esitniated. These 
events will be subtracted in determining the fi- 
nal branching fractions. Other surviving back- 
ground events mainly come from J /ip AE^tt"*", 
AS+TT", E+S-(1385), and E-S+(1385), and 
their charge conjugate channels, but they only 
give a flat contribution in the K^ signal region. 
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Figure 5. The missing mass spectrum of charged 
tracks in. J/ip —> nKgA + c.c. for events satisfying 
the requirements in Fig. 4 and Lxy{A) > 5 mm. 



The sum of these backgrounds, normalized by 
their branching fractions, is shown as the cross- 
hatched area in Fig. 6, and it is consistent with 
the background under the peak for data. 

Using a Gaussian to describe the Kg and a sec- 
ond order polynomial function to model the back- 
groimd shape, a fit to the m^+7r- distribution is 
performed, shown as the curve in Fig. 6. A to- 
tal of 1058±33 if Si events are obtained. No 
signal is observed in the m^+^r- invariant mass 
distribution for events which recoil against the A 
sideband region (1.140 < < 1.164 GeV/c^). 
The detection efficiency for the signal is 6.09%, 
which is determined from a uniform phase space 
MC simulation. The branching fraction is: 

B{J/il) nK%K + c.c.) 

^ Nobs - Nbg 

Nj/^ ■e-B{K^PK+)- B{K% ^ tt+tt") 

= (6.46 ±0.20) X 10"^, 

where Noi,g is the number of events observed 
(1058±33); N},g is the number of background 
events from J/ij: nKp:^ (42 ± 8) and pifgS" 
(12 ± 3); e is the detection efficiency; Njj^ is 




Figure 6. The 7r"'"7r~ invariant mass distribution 
for J/%1) nK^h. + c.c. candidates satisfying the 
A(A) mass selection requirements, Xic < ^' 
Lxy{A.) > 5 mm. The fit is also shown. The cross- 
hatched area is the sum of the backgrounds after 
normalization, as described in the text. 



the number of J/4> events; and B(A — > p7r+) and 

B(K'^ 7r+7T~) are the branching fractions of 
A pn'^ and Kg tt+tt" [6]. The error is 
statistical only. 

If we fit the charge conjugate channels sepa- 
rately, we obtain 502±22 events with an efficiency 
of 6.02% for J/ip nK%K, 560 ± 24 events, an 
efficiency of 6.16% for J/ip nK^K, and the 
following branching fractions: 

B{J/^ nK°A) = (3.09 ± 0.14) x 10"^, 
B{J/ip nJsfO A) = (3.39 ± 0.15) x 10"^, 

where the errors are statistical only. These results 
are consistent with each other in 1.5(T. 

In order to obtain a clean sample of J/tp 
nKgA and nKgA, we require events to satisfy 
the A(A) and Kg mass selection requirements, 
Xic < 5, and Lxy{A) > 5 mm, and also require 
the Kg decay length Lxy{Kg) > 5 mm to elimi- 
nate backgrounds without a Kg in the final state, 
such as J/V' — * AI]~TT"'". After this final selection, 
the background contribution is estimated to be 
less than 5%. Figure 7 shows the scatter plot of 



6 



for J/il) ^ nK^K candidate 
events for all but the A(A) and mass selec- 
tion requirements, where the boxes in the plot 
show the signal and sideband regions. The in- 
variant mass spectra of AiiT^, nK^^ and An(An), 
as well as the Dalitz plot for all selection require- 
ments are shown in Fig. 8. In the Aif^ invari- 
ant mass spectrum, an enhancement near h.K'g 
threshold is evident, as is found in the KK mass 
spectrum in J/i/; ^ pK~A [12]. If the enhance- 
ment is fitted with an acceptance weighted S-wave 
Breit-Wigner function and a function fbg{S) de- 
scribing the phase space "background" contribu- 
tion, the fit leads to M=1.648±0.006GeV/c2 and 
r = 61 ± 21MeV/c^, respectively. Here the errors 
are only statistical. The systematic uncertain- 
ties are not included since more accurate mea- 
surements of the mass and width should come 
from a full PWA involving interferences among 
A''* and A* states. The fitted mass and width 
are consistent with those obtained from a par- 
tial wave analysis of J/ip — > pK~A [12]. The 
X(2075) signal which was seen in the pA invariant 
mass spectrum in J/ip pK~A is not significant 
here. Using a Bayesian approach [13] and fixing 
the mass and width of X(2075) to 2075 MeV/c^ 
and 90 MeV/c^ respectively, the upper limit on 
the number of events observed N^^^ is 54 events 
at the 90% C.L. 

The N* state at around 1.9 GeV/c^ in the AK^ 
invariant mass spectrum and the A* states at 
around 1.5 and 1.7 GeV/c^ in the nK'g invariant 
mass spectrum are present. A larger data sample 
and a partial wave analysis are needed to analyze 
these A^* and A* states. 
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Figure 7. The scatter plot of m^-p versus m^+^- 
for J /ip nKgA candidates after all selection 
requirements except for the A and Kg mass re- 
quirements. The boxes in the plot show the signal 
and sideband regions. 



events. The statistical significance of the Kg is 
about 7.2(7. The 2 ± 1 background events from 
the A(A) sidebands and 2±1 background events 
from ip{2S) nif^E" are subtracted. A uniform 
phase space MC simulation determines the detec- 
tion efficiency to be 9.16%. The corresponding 
branching fraction is: 



B{^{2S) nK%A + c.c.) = (8.11 ± 1.14) x 10" 
Here the error is statistical only. 



3.2. Measurement of ^/j(2S') ^ nK^A + c.c. 

Using the same criteria as in Section 3.1, we 
select tp{2S) nKgA + c.c. events from the BE- 
SII sample of 14M ^{2S) events. The 7r"p and 
TT+TT^ invariant mass spectra, the scatter plot of 
rriTr-p versus 7717^+^-, and the missing mass spec- 
trum after the final selection are shown in Fig. 9. 
The A and K^ signals are obvious. 

Fitting the 7r"'"7r~ mass spectrum with a Gaus- 
sian for the K'g signal and a first order polynomial 
background, as shown in Fig. 10, yields 50±7 K'g 



3.3. Systematic Errors 

In this analysis, the systematic errors on the 
branching fractions mainly come from following 
sources. 

3.3.1. MDC tracking 

The MDC tracking efficiency has been mea- 
sured using channels like J /ip pw, AA, and 
ip{2S) ir+ir-J/iP, J/^ fi+^-. The MC sim- 
ulation agrees with data within 1 to 2% for each 
charged track [11]. Thus 8% is taken as the sys- 
tematic error coming from MDC tracking for the 
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Figure 8. The invariant mass spectra of (a) A.Kg, 
(b) nifg, and (c) nA, as well as (d) the Dalitz 
plot for candidate events after all selection crite- 
ria. The crosses show the sideband backgrounds. 



Figure 9. The invariant mass spectra of (a) tt~p 
and (b) tt+tt", (c) the scatter plot of M^-p ver- 
sus M^+^-, and (d) the missing mass spectrum 
for events satisfying the A (A) mass selection re- 
quirement, Xic < 5, and Lxy{A) > 5 mm. 
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Figure 10. The fit of Tr~^w~ invariant mass spec- 
trum for tp{2S) — » nKgA + c.c. for events satisfy- 
ing the requirements in Fig. 9d. 



4-prong events considered in this analysis. 

3.3.2. Kinematic fit 

The systematic error from the IC kinematic fit 

should be smaller than that from the 4C kine- 
matic fit, since there are fewer constraints. Var- 
ious studies show that the uncertainty of the 4C 
kinematic fit is around 4% [14] . Here we conserva- 
tively take 4% as the error from the IC kinematic 
fit. 

3.3.3. Particle identification 

In Ref. [11], the particle identification efficiency 
of TT, K, and p are analyzed in detail. Here, only 
one charged track is required to be identified as 
a p or p, and the systematic error from particle 
identification is less than 2%. 

3.3.4. A and Kg vertex finding 

In Ref. [15], J/ip ^ AA — > -!T~^Tr~pp is cho- 
sen as the reference channel to study the sys- 
tematic error of the A vertex finding algorithm, 
and 1.2% is determined as the systematic er- 
ror for one A vertex. For Kg, the efficiency 
of the secondary vertex finding is studied using 
J/tjj — > K*{892)K + c.c. events, and the system- 
atic error is about 4.1% [16]. 



3.3.5. MC model 

Different hadronization models for simulating 
the hadronic interactions give different detection 
efficiencies. Their differences are taken as sys- 
tematic errors. The systematic errors are 7.0% 
and 14.7% for J/iJ) nKgA and its conjugate 
channel, respectively, and 11.1% for tpi^S) — > 
TiKgA + c.c. The efficiency differences with or 
without considering the intermediate N* and A* 
states are also taken as the systematic errors. 
They are 5.3% and 4.5% for J/^ nK^A and 
J/ip — > nKgA, respectively. 

3.3.6. Background uncertainty 

The background uncertainties come from the 
uncertainties associated with the estimation of 
the sideband backgrounds, continuum events, and 
the events from other background channels, as 
well as the uncertainties of the background shape, 
different fit ranges, etc. Therefore, the statis- 
tical errors in the estimated background events, 
the largest difference in changing the background 
shape, and the difference of changing the fit 
ranges are taken as the systematic errors for the 
background uncertainty. 

3.3.7. Intermediate decay branching frac- 
tions 

The branching fractions of A ^ p7r~ and the 
Kg — > TT+TT" decays are taken from the PDG [6]. 
The errors on these branching fractions are taken 
as systematic errors in our measurements. 

3.3.8. Number of J/^; and ^(25*) events 
The total number of J/ip events is (57.70± 

2.62) X 10^, determined from inclusive 4-prong 
hadrons [17], and the total number of tp{2S) 
events N^(2S) is (14.0±0.6)xl0^, determined 
from inclusive hadronic events [18]. The uncer- 
tainty on the number of J/tp events, 4.7%, and 
the uncertainty on the number of tp{2S) events, 
4.0%, are also systematic errors. 

The above systematic errors are all listed in 
Table 1. The total systematic error is determined 
by adding all terms in quadrature. 
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4. Results 

The decays of J /ip and il){2S) to riK^gk+cc. are 
observed for the first time, and their branching 
fractions are: 

B( J/V' nK%K + c.c.) = 
(6.46_± 0.20 ± 1.07) X 10-^ 
B{J/^ nK^A) = (3.09 ± 0.14 ± 0.58) x 10-^ 
B{J/il) nK%k) = (3.39 ± 0.15 ± 0.48) x 10"^, 

B{ijj{2S) nKlK + c.c.) = 
(0.81 ±0.11 ±0.14) X lO"''. 

The ratio of the branching fractions of ^^{23) and 
J/tjj decaying to nKgA + c.c. is Qh = (12.6 ± 
3.5)% and obeys the 12% nilc well. 

There is no obvious enhancement near nA 
threshold. The upper limit on the branching frac- 
tion on the near-threshold enhancement X(2075) 
at nA threshold at the 90 % C.L. is: 

BiJ/tp K°sX) ■ B{X ^ nA + c.c.) 

~ Nj/^ • e • B(A ^ pn-) ■ B{Kl ^ TT+TT-) ■ (1 - Ssys) 

< 4.9 X 10"'"^ (90% C.L.), 

where iV^^ is 54 events; e=5.32% is the detection 
efficiency considering the angular distributions; 
N j/^ is the number of J/i]j events; B{A pir ) 
and B{Kg tt+tt") are the A pn~ and Kg 
7r+7r~ branching fractions, and 6sys is the sys- 
tematic error (17.3%). Taking into account the 
isospin factor, the branching fraction upper limit 
for B( J/V' K^X) ■ B{X nA + c.c.) is not in- 
consistent with that for B{J/ip KX) ■ B{X 
pA + c.c.) [4]. 
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Table 1 



Summary of the systematic errors. 



Sources 


Relative Error (%) 


decay modes 


J/i) nK^K 


J/V' ^ nK^A 


J/ip X{2075)K^ 


ip{2S) nJT^A + c.c. 


MDC tracking 


8.0 


8.0 


8.0 


8.0 


Kinematic fit 


4.0 


4.0 


4.0 


4.0 


Particle ID 


2.0 


2.0 


2.0 


2.0 


MC model 


15.6 


8.3 


9.4 


11.1 


Background uncertainty 


2.1 


5.1 


9.3 


7.1 


A and Kg reconstruction 


1.2 


1.2 


4.3 


1.2 


Intermediate decay branching fractions 


0.9 


0.9 


0.9 


0.9 


Numbers of J/tp and il}{2S) events 


4.7 


4.7 


4.7 


4.0 


Total systematic error 


18.9 


14.3 


17.3 


16.6 



